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Multishell helical gold nanowires (HGNs) suspended between semi-infinite electrodes are found to 
exhibit peculiar electron-conduction properties by first-principles calculations based on the density 
functional theory. Our results that the numbers of conduction channels in the HGNs and their 
conductances are smaller than those expected from a single-atom-row nanowire verify the recent 
experiment. In addition, we obtained a more striking result that in the cases of thin HGNs, distinct 
magnetic fields are induced by the electron current helically flowing around the shells. This finding 
indicates that the HGNs can be good candidates for nanometer-scale solenoids. 

PACS numbers: 73.63.Nm, 68.65.La, 73.21. Hb 
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As the miniaturization of electronic devices progresses, 
the fabrication of atomic-scale structures and the elucida- 
tion of their electronic properties have attracted consider- 
able attention. Within the last decade, a large number of 
experiments concerning electron-conduction properties of 
nanowires of a few nanometer lengths have been carried 
out using a scanning tunneling microscope or a mechani- 
cally controllable break junction Among the most 
exciting discoveries are the tip-suspended multishell he- 
lical gold nanowires (HGNs) of 1 nm diameter and 5 
nm length observed by Kondo and Takayanagi 0] using 
transmission electron microscopy (TEM). The HGNs are 
sandwiched between Au(llO) electrodes and composed 
of several atom rows twisting around their axis. In such 
minute structures, it has been well known that electron 
transport becomes ballistic and exciting quantum phe- 
nomena take place. For example, it was experimentally 
observed that as a gold point contact is stretched and 
thinned, the electronic conductance of the contact de- 
creases roughly in units of the quantized value, Gq — 
2e^//i, where e is the electron charge and h is Plank's 
constant When it comes to HGNs, since a few 
relationships between their geometries and conductances 
have been investigated using TEM images ||8|, there re- 
mains much to be learned about their electronic charac- 
teristics; in particular, the electron-conduction properties 
of these structures are important from both fundamental 
and practical points of view. Recently, Tosatti et al. 
explored the geometries and electronic states of infinitely 
periodic helical wires using conventional first-principles 
calculations with a plane-wave basis set. However, the 
conduction properties of nanowires intervening between 
electrodes cannot be determined from the results of in- 
finitely periodic wires because scatterings of electrons, 
which play crucial roles in the electron conduction, do 
not take place in the infinitely periodic systems. There- 
fore, in order to interpret the conduction properties pro- 
foundly by theoretical calculations, the models in which 



the nanowires are connected to semi-infinite electrodes 
are indispensable. 

In this study, we examined the electron-conduction 
properties of HGNs suspended between semi-infinite gold 
electrodes using first-principles calculations within the 
framework of the density functional theory. Our find- 
ings are that the numbers of conduction channels of the 
HGNs do not coincide with the numbers of atom rows 
in the HGNs, while the number of conduction channels 
in a single-atom-row nanowire is one. In addition, some 
channels of the HGNs are not fully open, whereas the 
channel transmission of the single-atom-row nanowire is 
close to one. As a consequence, the conductances of the 
HGNs are not quantized and become much smaller than 
those expected from a single- atom-row nanowire. More 
intrigingly, in some of the HGNs, electron currents ro- 
tate around the nanowire axis so as to induce a magnetic 
field. 

Our first-principles calculation method for electron- 
conduction properties is based on the real-space finite- 
difference approach [Tol |. which enables us to determine 
the self-consistent electronic ground state with a high de- 
gree of accuracy by the timesaving double-grid technique 
and the direct minimization of the energy functional 
|l2j |. Moreover, the real-space calculations eliminate the 
serious drawbacks of the conventional plane-wave ap- 
proach such as its inability to describe strictly nonpe- 
riodic systems. The local pseudopotential of a gold 6s 
electron constructed using the algorithm by Bachelet et 
al. 10 is employed to describe the ion-core potential, 
and exchange correlation effects are treated by the local 
density approximation We take a cutoff energy of 

12 Ry, which corresponds to a grid spacing of 0.91 a.u., 
and a higher cutoff energy of 110 Ry in the vicinity of 
nuclei with the augmentation of double-grid points. 

The computational models are as follows: The HGNs 
are sandwiched between semi-infinite gold electrodes im- 
itated by the structureless jellium rs=3.01 (see Fig. ^ 
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FIG. 1: (color) Isosurfaces of channel electron distributions of single-atom-row and 7-1 nanowires. (a) The first channel of the 
single-atom-row nanowire, (b) the first channel, (c) the third channel, and (d) the fourth channel of the 7-1 nanowire. The 
states incident from the left electrode are shown. The yellow and green spheres represent the atoms in the outer and inner 
shells, respectively. The rectangular parallelepipeds at both the edges of the nanowires are jellium electrodes. 



TABLE I: Numbers of atoms in the each shell of the nanowire. 





Innermost 


2nd 


3rd 


Single- row 


9 






7-1 


9 


63 




11-4 


24 


101 




13-6 


54 


121 




14-7-1 


9 


63 


133 


15-8-1 


9 


76 


139 



for example.). A unit cell of 40.1 a.u. in the x and y 
directions with periodic boundary conditions is adopted, 
where the x and y directions are perpendicular to the 
nanowire axis. The distance between jellium electrode 
surfaces is fixed at 51.2 a.u. The atomic geometries of 
the HGNs are determined usinsthe parameters evalu- 
ated by Kondo and Takayanagi Q , and their lengths are 
set at ^ 49.1 a.u. Table shows the number of atoms in 
the HGNs. The global wave functions for infinitely ex- 
tended states from one electrode side to the other are 
evaluated by the overbridging boundary-matching for- 
mula [15t llfi llTj l . The conductance of the nanowire sys- 
tem in the limits of zero temperature and zero bias is 
described by the Landauer-Biittiker formula G=tr(T'''T) 
Go 01 , where T is the transmission matrix. The eigen- 
channels are computed by diagonalizing the Hermitian 
matrix T^T figf. 



Table ITll lists the conductances and the channel trans- 
missions of the single-atom-row, 7-1, 11-4, 13-6, 14-7- 
1, and 15-8-1 nanowires suspended between gold elec- 
trodes. In the case of the single-atom-row nanowire, the 
number of conduction channels is one and the conduc- 
tance is quantized as 1 Gq, which are in good agree- 
ment with the results of previous theoretical calculations 
and experiments 4, 5, 6, 15, 20, 21] . In the case of a 
thicker HGN composed of several atom rows, the num- 
ber of conduction channels does not correspond to the 
number of atom rows, and some channels are not fully 
open. The reduction of the number of conduction chan- 
nels is mainly due to the geometry and electronic states 
of a nanowire rather than the interaction between the 
nanowire and electrodes; in fact, the number of channels 
in the 7-1 nanowire agrees with the number of energy 
bands crossing the Fermi level in the infinitely periodic 
wire 1^. On the other hand, the channel transmissions 
are closely related to the interfaces between the nanowire 
and electrodes; when we replaced the jellium electrodes 
with Au(llO) crystalfine electrodes, the conductance of 
the 7-1 nanowire becomes 4.8 Go, whereas the number 
of conduction channels remains six. Consequently, the 
conductance per atom row of the HGNs is much smaller 
than that of the single-atom-row nanowire. Compared 
with recent experimental results 1^, the calculated con- 
ductances agree with those of the experiments within the 
relative difference of 20 %. These results imply that the 
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TABLE 11: Conductances and channel transmissions at the Fermi level 





Single- row 


7-1 


11-4 


13-6 


14-7-1 


15-8-1 


Conductance (Go) 


0.96 


5.19 


9.08 


11.97 


13.82 


14.44 


1st ch. 


0.958 


0.997 


1.000 


0.996 


1.000 


1.000 


2nd ch. 


— 


0.995 


0.991 


0.992 


0.999 


0.996 


3rd ch. 


— 


0.970 


0.986 


0.982 


0.991 


0.991 


4th ch. 


— 


0.938 


0.890 


0.978 


0.984 


0.982 


5th ch. 


— 


0.653 


0.884 


0.962 


0.959 


0.975 


6th ch. 


— 


0.640 


0.874 


0.958 


0.947 


0.954 


7th ch. 


— 


— 


0.753 


0.940 


0.940 


0.942 


8th ch. 


— 


— 


0.748 


0.878 


0.927 


0.931 


9th ch. 


— 


— 


0.620 


0.872 


0.893 


0.899 


10th ch. 


— 


— 


0.517 


0.744 


0.888 


0.891 


nth ch. 






0.497 


0.735 


0.881 


0.841 


12th ch. 






0.318 


0.706 


0.867 


0.834 


13th ch. 








0.618 


0.770 


0.732 


14th ch. 








0.383 


0.700 


0.691 


15th ch. 








0.200 


0.542 


0.645 


16th ch. 








0.002 


0.527 


0.595 


17th ch. 












0.537 



electron-conduction channels of the HGNs are different 
from that of the single-atom-row nanowire; as an exam- 
ple, we show in Fig. ^ some of the channel electron dis- 
tributions of 7-1 nanowire, in which the states incident 
from the left electrode are depicted. For comparison, that 
of the single-atom-row nanowire is also described [Fig. 
(a)]. One can easily recognize that the conduction chan- 
nels do not completely form along the atom rows in the 
7-1 nanowire, while the conduction channel of the singlc- 
atom-row nanowire follows the atom row. Moreover, in 
Fig.n(d), the fourth channel of the 7-1 nanowire shapes 
on the outer shell and twists about the nanowire axis 
with a right-hand rotation, which is opposite to that of 
the atom rows. The helical channel current is mainly 
caused by the asymmetry of the atomic structure of the 
HGN. 

We next investigate the magnetic field induced by the 
electron current flowing through the HGNs. Since some 
channels twist around the nanowire axis, the helical elec- 
tronic current through the nanowires gives rise to the 
magnetic field inside the nanowire (see Fig. 12) ■ The direc- 
tion of the induced magnetic field approximately agrees 
with that expected from the channel electron distribu- 
tion. The magnetic flux density per bias voltage (lO^'^ 
T/mV) is of the same order as those obtained using the 
channel electron distribution and its transmission. It is 
worth noting that there is a considerable difference in the 
strength of the magnetic field between the thin and thick 
HGNs: the field in the thin nanowire is stronger than 
that in the thick one. In the case of the thick HGN, the 
axial symmetry of its atomic structure breaks very mod- 
erately, and then no distinct magnetic field is induced be- 
cause the electronic current does not twist largely around 



the nanowire axis [see Fig.|21(c)]. Okamoto et al. p3.l2^ 

attempted to explore the magnetic field using nanowires 
made of structureless solid jellium and reported that the 
strength of the magnetic field is very small. In the case 
of the jellium model, the shell structures of the HGNs 
are not taken into account, therefore the helical cur- 
rents flowing around the shells are not observed; the shell 
structures play important roles in the induction of mag- 
netic fields by the helical currents. 

To the best of our knowledge, the findings re- 
garding the electron-conduction behavior (conductance, 
conduction-electron distribution, and so on) of the HGNs 
and the magnetic field induced by the helical current 
are the first theoretical indications, and experimen- 
tal consequences should be interesting to pursue. Al- 
though whether the magnetic field would be measur- 
able in steady states similarly to the case of macroscopic 
solenoids remains an unsettled question, our findings in- 
dicate that the HGNs are good candidates for nanometer- 
scale solenoids. Thus, this situation will stimulate a new 
technology for electronic devices using atomic wires. The 
work in progress is a simulation of the formation of the 
magnetic field by the time-dependent density functional 
theory. 
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the 21st Century COE "Center for Atomistic Fabrica- 
tion Technology" and also by Young Scientists (B) (Grant 
No. 14750022) from the Ministry of Education, Cuffure, 
Sports, Science and Technology. The numerical calcula- 
tion was carried out by the computer facilities at the In- 
stitute for Solid State Physics at the University of Tokyo, 
and the Information Synergy Center at Tohoku Univer- 
sity 
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FIG. 2: (color) Induced magnetic flux density per bias volt- 
age, (a) 7-1 nanowire, (b) 11-4, and (c) 15-8-1 nanowire. The 
gray, yellow and green spheres represent the atoms in the in- 
nermost, second, and third shells, respectively. The broken 
lines are the edges of the jellium electrodes. 
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